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NEAR INFRARED-TRANSMISSION SPECTROSCOPY OF TONGUE TISSUE 



DESCRIPTION 

Background of the Invention 

This application relates to the non-invasive measurement and monitoring of 
chemical substances found in blood, and particularly to the non-invasive measurement and 
monitoring of blood glucose levels. 

Diabetes affects approximately 1 6 million people in the United States. People 
suffering from diabetes mellitus must monitor their blood glucose levels for insulin 
replacement therapy. At present, this generally requires an invasive procedure, performed 
several times daily, in which a sample of blood is taken for analysis. Non-invasive 
procedures have been proposed, including those found in US Patents Nos. 5,069,229 and 
5,459,317, and EP 0 548 41 8, which are incorporated herein by reference, but no in-home 
non-invasive glucose measurement device is as yet available. 

Among the proposals for non-invasive measurement of blood chemicals are 
those which depend on diffuse-reflectance of infrared radiation, for example from a finger, lip 
or forearm. In these techniques, light is delivered to the surface of the tissue, penetrates for . 
some distance and then is reflected back to be detected. The drawback of the diffuse 
reflectance measurements is that the depth to which the light penetrates is not deep enough to 
have enough interaction with glucose to provide consistent and quantitative detection of 
glucose. In addition, tissue areas conventionally used for diffuse reflectance measurements 
may have varying amounts of fat from individual to individual. Since fat interferes with 
detection of glucose, this variability can make calibration and quantification difficult. 

It is an object of the present invention to provide an improved method for the 
non-invasive measurement of glucose. 

Summary of The Invention 

The present invention overcomes the difficulties of prior art processes by 
using transmission spectroscopy, and by performing the measurements on the tongue of the 



v. 
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individual. The near infra red (NIR) first overtone spectral region between 7000 and 5000 
cm 0 (1 .43 - 2 urn) is well suited for non-invasive measurements. Because the tongue is one 
of the leanest accessible sights in the body, variations in readings due to body fat are 
minimal. Further, the tongue is well thermostated and vasculated which are important factors 
in selecting' a measurement site. 

Brief Description of the Drawings 

Fig. 1 shows SEP versus noise plots for model data sets. Vertical lines 
indicate noise levels of the webbing (solid) and tongue (dashed) for back to back spectra at 
different measurement wavelengths. 

Figs. 2A-F show absorbance spectra for near-infrared measurements at 

different measurement sites. 

Fig. 3 shows single beam spectra of human webbing (dashed) and tongue 

(solid) with a glucose absprbance spectrum (dotted) overlaid. 

Figs. 4 A and B show 100% lines of back-to-back psectra of human webbing 
and tongue with an absorbance spectrum of 1M glucose overlaid. 

Figs. 5A-E show glucose concentrations versus sample number for the 
calibration samples for five individuals. 

Figs. 6 A-D show SEP versus RMA noise calculated over 100 cm' 1 spectral 
ranges for model data sets. Circles and triangles are 5.6 and 6.3 mm aqueous thickness 
without scatter, respectively. Squares are prediction errors from the model data set with 
scattering particles added. The vertical lines are noise levels of subject D who had an 
aqueous layer thickness of 5.5 mm. 

Figs. 7A-E are concentration correlation plots for the best PLS models for the 

five volunteers predicting the blind samples. 

Figs. 8A-E are concentration plots for the five volunteers, using the first 
sample as the calibration and the last samples as the prediction. 
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Detailed Description of the Invention 

In accordance with the present invention, non-invasive monitoring of glucose 
in a human patient, is carried out by measuring the near infrared absorbance of glucose using 
transmission spectroscopy on the tongue of the human patient. This selection of 
measurement site provides better correlation with invasive measurements than non-invasive 
measurements performed at other sites, including the index finger and webbing of the hand. 

In selecting a site for non-invasive glucose determination, several factors need 
to be considered. First, a useful technique needs to produce results which are reliable 
predictors of the actual glucose levels. Second, it would be desirable to have a calibration 
model which would work in substantially all patients. Patient-to-patient variation in the 
physical size of the tissue being measured and in the levels of interferents such as fats, which 
can act as a major interferent in the measurement of glucose in the overtone region, can make 
achieving these goals a complicated process. Our evaluation of several alternative 
measurement sites, however, indicates that the tongue is a highly suitable location for non- 
invasive glucose measurements. 

To evaluate various sites for their suitability, single beam transmission spectra 
were collected of human webbing between the thumb and forefinger, tongue, nasal septum, 
cheek, upper lip and lower lip. As described below, the hand webbing had the highest levels 
of fat in the individuals tested, thus producing the worst results, while the tongue had the 
lowest fat levels and produced the most accurate and reproducible results. 
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EXAMPLE 1 
Evaluation of Alternative Measurem ent Sites 
Spectra were collected on a modified Midac spectrometer for the near infrared 
range. Modification include laser detectors, a fan designed to control temperature in the 
spectrometer, a 250-Watt source operating at about 1 10 watts and a one-millimeter diameter 
InGaAs detector with a 1 .9 urn cutoff (Epitaxx). An H-band astronomical filter was used to 
limit the source bandpass to utilize the full dynamic range of the detector without saturation 
and minimize sample heating. A one-inch diameter, 25 mm focal length convex lens was 
used to focus the source intensity on the sample or fiber optic bundle. 

Preliminary single beam spectra were collected using a fiber optic bundle (Dolan 
Jenner Industries) connected to the thickness controlling window to direct the incident 
radiation to the sampling site. The detector and mount was freed from the spectrometer 
for case of sampling. This arrangement allowed spectra to be collected from hard to 
reach body sites. Single beam transmission spectra were collected of human webbing 
. between the thumb and forefinger, tongue, nasal septum, cheek, upper lip, and lower lip. 
A background single beam spectrum was collected without a sample in the optical path 
using a combination of a neutral density filter to limit the incident source intensity and a 
low detector preamplifier gain setting to avoid detector and electronic saturation. The air 
spectrum was corrected to match the higher gain and incident radiation intensity of the 
tissue spectra before converting it to absorbancc. 

Spectra with lower noise levels were collected without the fiber optic bundle to 
maximize incident source intensity. The detector and mount were fixed in the 
spectrometer in a vertical position for the tongue spectra and in a horizontal position for 
the webbing spectra. A mirror was used to direct the incident radiation to the detector for 
the tongue spectra. 
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Absorption spectra of the webbing of 19 volunteers were collected. A regression 
analysis was performed to compute the thickness of water and fat from the spectra using 
beef fat and water as the standards. The average thickness of water was 4.67 mm and 
beef fat was 2.4 mm for a webbing thickness of 6 mm. Table 1 lists RMS noise levels 
for various frequency ranges of ten back-to-back webbing spectra of 128 scans each. 
Figure 1 shows plots of SEP versus RMS noise levels for various frequency ranges for 
tnodcl spectra consisting of a constant fat thickness of 1 .6 mm and variable water 
thi ckne S sesof5.6and6.2m m .From Figure 1 and Table 1 , a very high SEP would 
be expected if the webbing were used (greater than 3 mM). The RMS noise levels are 
low in the spectral region which is dominated by water absorbance (>6000 cm"'). The 

of magnitude higher than the model spectra. Between measurement variations in the 
amount of fat compressed within the light path have not been successfully controlled. 
Another body site for the noninvasive measurement is needed because of the strong 
overlap between the absorbance bands of fat and glucose. 



Table 1 RMS noise values calculated from webbing and tongue 100% 



Spectral Range (cm 1 ) 
Sample 


6400 
-6300 


6300 
-6200 


6200 
-6100 


6100 
-6000 


6000 
-5900 


5900 
-5800 


6 mm webbing 
5.5 mm tongue 


54.7 
96.7 


25.1 
47.7 


18.0 
25.4 


36.1 
18.0 


271.1 
56.6 


605.8 
68.! 
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The composition and physical thickness of several bodily sites were evaluated 
with the goal of finding a suitable site for a noninvasive measurement of glucose. 
Absorption spectra are listed in Kg"" *A through 2F. It is apparent from the 
spectra that the upper and lower lips, tongue, cheek and nasal septum have less fat than 
the webbing. A regression analysis was performed to quantify the amount of fat and 
aqueous thickness. Table 2 lists the calculated aqueous and fat thicknesses. The 
tongue was found to be most desirable because it has the least amount of fat and there are 
relatively few variations between people. Major disadvantages of the tongue are the 
possibility of spreading disease while sampling many people and the discomfort of 
collecting spectra. 

Table 2 Calculated aqueous and fat thicknesses obtained from regression analysis. 



Sample 


Thickness 
(mm) 


Aqueous thickness 
(mm) 


Fat thickness 
(mm) 


Scatter 


Webbing 


5.75 


4.91 


2.62 


1.56 


Upper lip 


6 


5;91 


0.48 


1.72 


Lower lip 


6 


5.95 


0.65 


1.84 


Tongue 


6 


5.61 


0.26 


1.47 


Nasal septum 


6 


5.86 


0.51 


1.81 


Check 


6 


5.61 


0.37 


1.61 
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Spectra were collected through the tongue of various people. The average 
thickness of fat and water was 0.2 and 5.9 mm, respectively, for the 10 volunteers who 
were sampled. Figure 3 shows single beam spectra of the human tongue and webbing 
collected under the same instrumental configuration with a superimposed absorbance 
spectrum of 1 .0 M glucose. Although the spectrometer alignment was not held constant, 
the spectrum of the tongue is less absorbing and has less fat. Figure 4 shows 100% 
lines of back-to-back spectra of the human webbing and also from human tongue. Urge 
variations in the fat bands which overlap with the glucose bands at 5920 and 5750 cm' in 
the webbing spectra makes it obvious that the tongue is a superior measurement site. 

Table 1 also shows noise levels of tongue spectra calculated in the hospital with 
multiple interfaces. From these noise levels and Figure 2 one may estimate a SEP of 
about 2 mM on a one subject study. 

It is clear that the tongue is a superior measurement site for a noninvasive 
measurement compared to the webbing. The convenience of collecting spectra through 
the webbing of the hand is traded for lower variations in the amount of fat within the light 
path. The interface between the spectrometer and the volunteer will need to be sterilized 
between collecting spectra from different subjects. If the same spectral quality were 
obtained in a noninvasive data collection as the test spectra from the tongue a SEP of 
about 2 could be obtained. 
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EXAMPLE2 

Bgtgctjon of BlooH glucose. Levels T icing Iransmjssion ^prr-tr nsropy of The Tongue 

Spectra were collected over the course of 39 days on five type 1 diabetic 
volunteers, 3 female and 2 male. The experimental procedure was approved by a human 
subject protocol board attheUnivcrsity of Iowa Hospitals and Clinic, Two spectra and 
one reference blood glucose concentration were collected for each volunteer's visit which 
numbered greater than 100. The number of visits per day were limited to less than five 
per volunteer. 

Spectra were collected using a Midac M-Series specrrometer modified for the 
m. The spec.rpme.er modifications arc detailed in Chapter 8. The source was held * 
8 1 .0 watts by supplying constant voltage. An H-band as.ronomical filter was used to 
Mtthcbandpass of thesourcc intensity The l-mm diameter 1.9 pre cutoff InGaAs 
dereetor (Bpitaxx) was held at 1 5 K by supplying constnnt vo.rage to a thermal electric 
code, The resistor in .he feedback loop of the firs, stage of me doctor preamplifier 

UX mterferograms <8K, worn col.ec.cd. Fourier .reformed, and the resulting 
single beam spectra with a da.a spacing of 4 cm- were saved over the spcotra, range 
borween 7000 and 5000 cm" 1 . The .ctnperatnte of the interface was held a. 38.5 «C using 
i DC hcting element and a temperature con.ro. cireui,. The .emperauare was monitor 
at the following oites during the course of each data colleofiont room, interfaoe. and 
inside the spectrometer. 
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Five interchangeable stainless steel interfaces were machined. These consisted of 
a cylindrical sheath to cover. the detector and an adjustable lever which controls the 
sample thickness. Each piece had a hole machined through with a 1 1.5-mm diameter 
circular sapphire window to allow the source radiation to pass through the sample to be 
detected. The thickness of the sample was confined to 5.45 mm. The interchangeable 
interfaces were swabbed with 95% ethanol and autoclaved at 130 °C for 5 minutes 
between volunteer visits. 

Reference Blood Glucose Values 
Capillary blood from a finger stick was drawn into a 100 |iL volume heparinized 
tube. These samples were analyzed immediately with a YSI 2300 Stat Plus glucose 
analyzer (Yellow Springs Instruments) with both probes measuring glucose. These 

reference values were validated with a simultaneous measurement using a One-touch test 
strip capillary blood glucose analyzer. The average of two YSI glucose values was used 
to build PLS models. Every fifth blood sample was considered a 'blind' and was not used 
. to optimize the PLS models. 

Data Analysis 

Separate PLS models were optimized using single beam spectra from each 
individual. The first 189 spectra were not included because the source voltage varied and 
the volunteers were learning how to perform the measurement in this time period. 
Spectra were also removed if their single beam intensity was very low. A PC A score 
analysis was performed but no spectral outliers were removed. The PLS calibration 
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model parameters were optimized. Three unique rearrangements of calibration and 
monitoring sets were used to optimize spectral range and number of PLS factors. The 
parameters which gave the lowest average SEM for all three rearrangements of calibration 
and monitoring were selected. This procedure reduces the chance of over modeling and 
modeling information specific to a single small monitoring set . 



Concentration Correlations 
By limiting the number of samples that were collected each day, the correlation 
between glucose and time was minimized. . Figure 5. shows plots of glucose versus 
sample number for the calibration samples for the five subjects. They look like scatter 
plots. Table 3 shows the correlation coefficient (r 2 ) of the time profiles. The 
correlation coefficients are all below 0.0 14. It is not likely that any temperature variation 
or spectrometer drift correlated with these randomized concentrations. In this way, 
spectral variations that the PLS models are built upon should only be analyte specific 
changes. 



Table 3 Description of glucose concentrations. 



subject number high low mean 
(mM) (mM) (mM) 



A 
B 
C 
D 
E 



173 
178 
181 
207 
169 



21.3 
28.3 
19.1 



2.94 
2.54 
2.83 



26.85 3.20 
16.75 3.79 



9.91 

12.64 

9.05 

11.88 

8.57 



•standard deviation 
(mM) 



4.53 
6.67 
4.06 
5:09 
2.83 



Regression coefficcnts of time plots 



Bo 



Bi 



8.904 0.01082 

12.13 0.006442 

12.62 -2.203e-4 

8.625 5.032c-3 
8.793 -8.642c-4 



0.01377 

2.465c-3 

6.085c-6 

4.245c-3 

2.020c-4 
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S-.gnif.cant glucose concentration variation is seen in all subjects except volunteer " 
E who has very tight control over her giycemia. Table 3 Iists the high, low, mean, and 
standard deviation of each subject's sample concentrations as well as the total number of 
spectra. The subjects were not persuaded to vary their daily diabetes treatments. The 
gl ucose concentrations are be.ieved to be representative of the blood glucose variation for 
these individuals. Significant spectra, and concentration changes are needed to build 
meaningful PLS calibration models. The limited concentration range for volunteer E 
might make calibration difficult. 

Spectral Quality 

The spectral quality is fairly high. Table 4 lists RMS noise levels calculated 
over 100 crn 1 intervals, the standard deviation of the ratio of the intensity at 5751 to 6994 
cm' 1 , standard deviation of the single beam intensity, and signal to noise ratio. The RMS 
noise levels were calculated using a second order polynomial fit for 100% lines from 
back-to-back single beam spectra. The noise levels between 6100 and 6000 cm* 1 are 
fairly low but the noise levels between 6000 and 5900 cm" 1 are not as low as those 
collected on the author in the pre-cxperiment trial. The ratio of the intensity at 5751 to 
• 6994 cm"' is the intensity at the fat band divided by the intensity at the peak of the single 
beam, respectively. The signal to noise (SNR) is the peak single beam intensity divided 
by the RMS noise of the single beam calculated between 7000 and 6800 cm -1 . The 
SD/mean intensity is the standard deviation of the peak single beam intensities divided by 
the mean peak intensity multiplied by 100%. 
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Comparison to Model 
The in vitro model was used to estimate ideal prediction errors for the noninvasive 
data set. The aqueous layer thickness was estimated to be about 5.5 mm. By comparing 
the noninvasive noise levels to Figure 6 , it was estimated that an SEP of about 2.3 mM 
would be the best that could be expected for the noninvasive data if there were no 
interferences and no variation in the amount of fat in the optical path. The vertical lines 
are the noise levels from subject D. The SEP estimate for the 6000-5900 cm"' plot is 
found by drawing a line between the models with and without scattering particles. 



Table 4 Description of noninvasive spectral quality. 



Subject 



A 
B 
C 
D 
E 



6100-6000 
(UA.U.) 



26.9 
28.0 
26.6 
24.5 
31.1 



6000-5900 
(UA.U.) 



95.1 
88.3 
89.8 
82.0 
99.1 



SD Intensity , SD/mcan 
(575 1/5994 cm" 1 ) inlcnsity(%) 



0.007191 
0.009485 
0.005072 
0.006365 
0.011918 



7.71 

8.08 

8.34 

9.64 

9.62 



SNR 
(xlO 4 ) 



4.2 
3.5 
4.2 
3.9 
3.7 
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Predicting Blind Samples 
The calibration data set was used to predict the blind samples using the optimal 
PLS parameters. Figure 7 A-E are concentration correlation plots for the best 

PLS models for the volunteers predicting the blind samples. The circled points are 
spectra in the blind that are not within the scatter of the calibration data. The circled 
point for volunteer D at about 1 1 mM was had the highest normalized intensity at 5943 
cm" 1 which is where protein absorbs. The circled point for volunteer C at 19 mM is the 
highest concentration sample in the data set. There seems to be a bias of under predicting 
high concentration samples. This bias was also observed in monitoring sets. The best 
SEP's are 2.70 and 2.66 for volunteers D and C with the high concentration samples 
removed, respectively. 



Table 5 °P tiraal PLS parameters and prediction errors from predicting blind 
samples. 



SEP, includes all data in blind 
SEPb excludes circled data 
* SEP is greater than SDC 



Who 


• Range 
(cm*) 


Factors 


SEC 
(mM) 


SEP, 
(mM) 


SEP b 
(mM) 


SDP 
(mM) 


A 


6510-5600 


II 


3.51 


3.48 




4.41 


B 


6310-5550 


11 


4.29 


7.56* 


.5.50 


6.75 


C 


6600-5650 


11 


3.04 


3.49 


2.66 


3.73 


D 


6550-5620' 


16 


2.52 


3.96 


2.70 


4.69 


E 


6600-5650 


10 


2.62 


3.35* 




2.74 
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,„ Figure 7B for subject B one point is not shown. The mode, predicted a blind 
rop ,e ,4S.9 mM when the actua, con— was ,6.5 mM. The SEP without this 

model for subjecl E is predicting the mean concentration. The 
correlation plot (Figure >7E). Even 



point is 5.50 mM. The i 
slope is actually negative for this concentration 
lh „„gh the model has no prediction, 88% of the prediction points are in the A and B 
rcg io„s due ,0 the limited concentration range of this subject. • Table 6 lias the SEC, 
SEP, regression nna IyS is and percentage of predicted samples in each of the regions on 
the Clarke error grid. 



,«ri Clarke error grid analysis of predicting blind 
Table 6 Prediction errors and Clai Kc ei roi y 
1 " ' ~ samples. 



subject 


SEC 
(mM) 


SEP 
(mM) 


SDP 
(mM) 


A 


3.51 


3.48 


4.4! 


B 


4.29 


7.56* 


6.75 


C 


3.04 


3.49 


3.73 


D 


2.52 


3.96 


4.69 


E 


2.62 


3.35* 


2.74 



Bn 



B, 



4.031 0.5466 0.4333 54.1 

,0.01 0.3393 0.1198 28.9 

5.423 0.3710 0.2274 43.2 

5.992 0.4543 0.3136 53.7 



B 


C 


D 


E 


37.8 


8.1 


0 


0 


52.6 


13.2 


5.3 


0 


51.4. 


0 


5.4 


0 


34.1 


4.9 


7.3 


0 


44.1 


0 


11.8 


0 
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Predicting Ust Part of Data 
I„ a real world situation for in-home or clinical devices there will be an initialdata 
collection used to build the calibration set. The model must be robust so it can accurately 
predict the concentrations of spectra following the initial calibration. Predicting the blind 
samples which were every fifth sample throughout the entire data collection is unrealistic 
in that there is a calibration sample before and after each blind sample. Spectrometer or 
temperature drifts could be more easily accounted for in this situation. To approximate a 
real world calibration, spectra from the beginning of the data collection were used as the 
calibration and spectra from the end of the data collection were predicted. The same PLS 

. Fieures 8 A through E are 
parameters used to predict the blind samples were used. * 8 

concentration correlation plots of the five volunteers using the first samples as the 
calibration and the last samples as the prediction. Table 7 shows details of the data 
sets and model performance. 



Table 7 



Prediction errors 
data. 



and Clarke error grid analysis of predicting last part of 



subject 


SEC 


SEP 




(mM) 


(mM) 


A 


3.27 


4.42 


B 


4.77 


6.08 


C 


2.96 


3.86 


D 


2.73 


3.35 


E 


2.71 


2.72 



Regression Analysis of % of prediction in Clarke error grid Regions 



Prediction Data 

n n- 



4.11 8.027 0.3358 0.3067 40.7 
5.44 4.119 0.6930 0.6303 54.2 35.4 . 2.1 



B 


C 


D 


E 


55.8 


5.8 


11.5 


0 


52.9 


20.6 


2.9 


5.9 


50 


5.6 


1.9 


1.9 


35.4 


. 2.1 


10.4 


0 


56.2 


0 


0 


0 
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In the noninvasive spectra there are two sources of RMS noise on 100% lines. 
The first is spectral noise from the detector and other electronic and optical components 
of the spectrometer. This noise is related to the number of photons at a given frequency 
which are detected by the detector. The in vitro, model of a aqueous and fat layer of 
constant thickness produces this type of noise. The other type of noise is variation of 
strongly absorbing components such as fat and muscle. Although the noninvasive spectra 
have low noise in the spectral region dominated by water absorbance (greater than 6000 
cm-'), there are variations in the fat and muscle bands (below 6000 cm" 1 ) which are 
believed to be causing our high prediction errors. Subject E who had the worst PLS 
mo dc. had the highest variation in the fat band. This variation could be related to body 
composition or movement during the data collection. 

The concentration correlation plot for. subject E illustrates a situation where all 
predictions are in the A and B regions on the Clarice error grid but the model is clearly 
void of any analytical utility. This shows that the Clarke error grid alone is insufficient to 
judge analytical performance of PLS models. 

Although, it is debatable whether glucose is being modeled for subjects A, B. and 
C. There is very strong evidence that glucose is being modeled for subject D. There is a 
bias of under predicting high concentration spectra. If this were due to a lack of 
analytical information the lower concentrations would be over predicted as with subject 
E. This is not the case. The predictions at low concentrations actually follow the unity 
line. It has been reported that there is a time delay in the increase of interstitial glucose 
concentration when blood glucose concentration is rising. When light passes through 
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the tongue it is-belicvcd that the majority of the photons will internet with interstitial fluid 
not blood! The mismatch between our predicted glucose values and reference blood 
glucose values at high concentrations accounts for scatter at high concentrations. There 
are conflicting claims in the literature about the exact relationship between blood and 
interstitial glucose concentrations. If the lower interstitial glucose does not track blood 
glucose at high concentrations, it would be further evidence that we are actually 
measuring glucose. This has never been observed before and is further evidence that we 
are measuring glucose absorbancc information. 

For the first time, a human noninvasive blood glucose measurement has been 
demonstrated. The calibration was such that glucose concentration was independent of 
sample order which eliminates the possibility of falsely low prediction errors due to 
chance correlations between glucose concentration and temperature or spectrometer drift. 
A real world calibration where data collected initially was used to predict data collected at 
a later time yielded similar prediction errors as predicting the blind samples which were 
every Fifth sample throughout the data set. 

Future «ements, an improved design for the interface may improve 
spectral reproducibility, thereby, reducing prediction errors. Also, the source power could 
be increased by at least 50% which would yield a higher signal to noise ratio. Although, 
spectra, quality was probably compromised, the two-month data collection better 
approximates real world situations compared to calibrations over the course of one week. 
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CLAIMS 

1 . A method for non-invasive monitoring of glucose in a human patient, 
comprising measuring the near infrared absorbance of glucose using transmission 
spectroscopy, wherein the measurement is performed on the tongue of the human patient 

2. The method of claim 1 , wherein the near infrared absrobance is 
measured in the first overtone spectral region between 7000 and 5000 cm"'. 
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